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Abstract —In this research, total acceleration equation is 
formulated where there is time scale coefficient at its time 
differential term. The formulation was done based on 
Courant Number equation and by using Taylor series. 
Then this total acceleration is applied to kinematic free 
surface boundary condition and Euler momentum 
equations. Potential velocity and water surface equations 
of linear water wave theory as well as wave number 
conservation equation were substituted to momentum and 
kinematic free surface boundary condition equations 
produced dispersion equation with wave amplitude as its 
variable and which fits with wave number conservation 
equation. Wave number consen’ation equation is an 
equation that regulates changes in wavelength as a result 
ofwaterdepth changes. This equation was extracted from 
potential velocity equation. 

Keywords—Courant Number, Taylor Series, Total 
Derivative, Wave Number Conservation Equation. 

I. INTRODUCTION 

This research aims at finding a dispersion equation with a 
wavelength that fits with what exists in the nature. 
Dispersion equation of linear wave theory (Dean, 1991), 
was formulated using kinematic free surface boundary 
condition and Bernoulli equations where this Bernoulli 
equation is formulated from Euler momentum equation. 
Both at the kinematic free surface boundary condition and 
Euler momentum equations there are total change term of 
spatial and time function. 

In the formulation of total acceleration equation, it is 

dx 

defined that at limdx and dt close to zero, u = — where 

dt 

u is a velocity of material movement. Courant number in 
fluid mechanics (1928) stated that in order to be defined 

dx 

that u = —, there are certain criteria for the size of space 
length (dx) and time step (dt). In general, it can be stated 
that in order to be defined that at limdx and dt close to 

zero, u = —, there are certain conditions. 

dt 

The accuracy of Taylor series is determined not only by 
the number of its terms but also by its interval size. 
Meanwhile, Taylor series is often used only up to the first 


derivative where in this case, tmncation error can only be 
reduced by reducing the size of its interval. Based on 
Courant criteria it is assumed that Taylor series for a 
spatial and time function must contain a coefficient at the 
time interval. Then minimizing truncation error at Taylor 
series obtains time coefficient value and interval size that 
produces the level of accuracy that fits with what has been 
demanded. 

H. THE FORMULATION OF TOTAL 
ACCELERATION EQUATION WITH 
COEFFICIENT 

In this section total acceleration equation will be 
formulated where there is a coefficient at the time 
differential term. 

2.1. Base of the Theory 

a. Courant Number 

Courant (1928) introduced Courant Number which is a 
criteria relation between length interval (bx)with time 
step (St) to conduct numerical analysis at the fluid flow, 
i.e. 

C = — < c (1) 

where itis a velocity, <Stis time step and bxis length 

Sx 

interval, C max = 1. If at (1) it = —is defined, hence C = 1 

which does not meet the Courant Number criteria. 
Sx 

However, if u = —is defined where y is a positive 

number greater than 1, then it will meet (1). From this 
equation, a conclusion can be made that there is a 
coefficient y at time step St to define a velocity. This 
coefficient can be stated as a time scale coefficient. 

b. Taylor Series Review 

Taylor series is often used only up to the first derivative 
or with an accuracy of 0(5 1 )at numerical analysis as 
well as the formulation of a conservation law. Total 
acceleration equation, at Euler momentum equation at 
fluid flow, is also often formulated using Taylor series 
0 (5 1 ) Using Taylor series up to the first derivative, the 
accuracy is depended only on the interval measurement. 
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Therefore, interval size that produces a good accuracy 
should be determined. 

For a continuous function f = fix, t), where xis 
horizontal axis and tis time, Taylor series approach 
OCd 1 ) (Thomas (1996)), 


fix + Sx, t + St) = fix, t) + 8x — +8t — 
J J ax at 

Using time scale coefficient y, 
f(x + Sx, t + y8t) = fix, t) + Sx + ySt 
.( 2 ) 

At (2) there is a tmncation error, 

„ Sx 2 a 2 / , y 2 st 2 a 2 / , r ^ a 2 / , 

R = --H-- + Sx ySt -h ••• 


a / 


2 ax z 


at z 


3x31 


R can be ignored if 


~TT 

c ax 4 


< £ 


..(3) 


where Eis a very small number. (3) can be achieved by a 
small size of dxandbtand with a time scale value of y. 
The size ofbxandbtand the value of y where Taylor series 
can be used only up to the first derivative can be 
determined with (3). 


2.2. Determining 8x,St andy. 

With an assumption that there are the values of 
5x,5tandy, where the number of the third derivative term 
with higher derivatives is much smaller than the number 
of the second derivative term, then as R, only the second 
derivative that can be used, and (3) becomes 


Sx 2 a 2 f v 2 st 2 a 2 f 


2 ax 2 


2 at 2 

AT. 


ySxSt-t 


Sx ^ + y s ^ 


< £ 


■ ■(4) 


Wave length is L = CTwhereC is wave celereties or wave 
velocity and T is wave period. From this wave length 

Ott 

equation, Sx = CSt 

- is obtained, then, 

K 

Sx = — .(5) 


defined. With L = —where k is 

k 


wave number, relation C 

o 8 t 
k 

Substitute (5) to (4) and the upper and lower part of the 
equation are divided by St 
i(a \ 2 ... a 2 f. to„.\ a 2 f . v 2 sta 2 f 

zW 5t ^ + n? 5t W +1 T-^ 


Vfc'a 


fc' ax +Y at 

Furthermore a sinusoidal function is reviewed with the 
form/ (x, t) = coskxcosot. This equation is a water 
wave surface elevation equation of the linear wave theory 
(Dean, 1991). Substitute fix,t) = coskxcosot to (6) and 
perform it at the condition coskx = sinkx = cosot = 
sinot, will produce 


< £ ....(6) 


—-+ y-^—\(rSt 
I 2 ' 2 I 


(l+y) 


<£ ....(7) 


At (7), the lower part of equation can be taken out from 
the absolute operation I I, since it always has positive 
value. At (7) there are two unknowns, i.e. yandbt. 


Another equation is needed, therefore Taylor series is 
performed with a change in t time only. 


fix,t + 5t) = fix.t) + St —+-?-{+■■■ ( 8) 

J J dt 2 dt 2 


In this case the time scale coefficient y was not performed 
since the one to be reviewed is only the change in 
function against time t. In order for (8) to be able to be 


performed only up to the first derivative, then 


St 2 d 2 f 
2 dt 2 




< £ ....(9) 


Substitute/Gc, t ) = coskxcosotto (9) and the equation is 
performed at coskx = sinkx = cosot = sinot, hence 

8t ma x< 7 ....(10) 

btat (10) is the value of St max since it is determined only 
based on the function of time without the interaction with 
the changes toward space. Withbf = St max , the value of y 
can be calculated at (7). 


Table.1: The Value ofyandSt max for fix, t) = 
coskxcosot. 



Y 

^max 

Wave period T: 1 sec. 

R 2 

3 

0,02228 

R 3 

3,14301 

0,02243 

R 4 

3,15978 

0,02214 

Wave period T: 8 sec. 

R 2 

3 

0,02547 

R 3 

3,14301 

0,02564 

R 4 

3,15978 

0,0253 

Wave period T: 9 sec. 

R 2 

3 

0,02865 

R 3 

3,14301 

0,02884 

R 4 

3,15978 

0,02846 

Wave period T: 10 sec. 

R 2 

3 

0,03183 

R 3 

3,14301 

0,03205 

R 4 

3,15978 

0,03162 


Table (1) presents the result of the calculation of values 
yand8t max for various wave periods and various level of 
accuracy R, where)? 2 shows that I?is calculated only with 
the second derivative only, R 3 shows thatf?is calculated up 
to the third derivative and R A , Ris calculated up to the 
fourth derivative. It can be seen that for the same period, 
the higher the accuracy of f?the bigger the value of y but 
with small change. Whereas at the similar accuracy level 
of R, for different wave period, the value of y is the same. 
The uses of accuracy up to R 5 does not change the value 
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of yconsidering the terms of and—at!? 5 is a very small 

‘ D 120 120 J 


number close to zero. 


The function fix, t) = coskxcosat has been used to 
calculate the value y. If the form of function fix,t) = 
sinkx sinat is used, where this equation is the change of 
water particle velocity at space and time for standing 
wave, then relation equation between yanddtis obtained 
which is similar to (7) for accuracy R 2 , i.e. 

——+ y-*—\(T8t 

d +/ ) < £ . 

However, accuracies/? 3 a nd/? 4 have different shapes and 
produce different value y, although with a not too big 
different, as shown on Table (2). 


Table. 2: The Value ofyandSt max forfix, t) = sinkxsinijt 



Y 

max 

(sec.) 

Wave period T = 7 sec. 

R 2 

3 

0,02228 

R 3 

2,85619 

0,02243 

R 4 

2,87499 

0,02214 

Wave period T = 8 sec 

R 2 

3 

0,02547 

R 3 

2,85619 

0,02564 

R 4 

2,87499 

0,0253 

Wave period T = 9 sec 

R 2 

3 

0,02865 

R 3 

2,85619 

0,02884 

R 4 

2,87499 

0,02846 

Wave period T = 9 sec 

R 2 

3 

0,03183 

R 3 

2,85619 

0,03205 

R 4 

2,87499 

0,03162 


Table (2) shows that at fix, t ) = sinkxsinat, the value of 
y is fluctuating against the level of accuracy /?but with 
relatively small fluctuation. From the two analysis of 
coefficient y for the two shapes of the function, the 
hydrodynamic analysis for water wave can use the value 
of y = 3. For numerical analysis where discretization of 
space and time is needed, than the space length size Sx = 

-with tune stepot =-. 


2.3. Total Acceleration with coefficient 
As has been shown that by performing coefficient on time 
differential term, Taylor series can be performed up to the 
first derivative, i.e. 


fix + Sx, t + St) = fix, t) + ySt — + Sx — 

’ ’ at ax 

The first temi of the right side of the equation is moved to 

the left and the equation is divided by St 

fix + Sx, t + St) — fix, t) 3 / Sx 3/ 

St ^ 3 1 St dx 

With the presence of time coefficient yat time differential 
term, it can be defined that at limSx, /(/approaches zero 
Sx 

can be defined that— = u. Therefore, the total 
st 


acceleration equation is 
of 


a/ a/ 

— — y -1 - u — 

dt a t 3x 


( 12 ) 


III. EQUATIONS FROM VELOCITY 
POTENTIAL 

This part has been written by Hutahaean (2019), however; 
considering that equations in this part are very important 
for this research, it will be rewritten. 


3.1. Velocity Potential Equation 

Velocity potential from linear wave theory which is the 
product of Laplace equation operation (Dean, 1991) is 

<t>ix,z,t ) = Gcoskxcoshkih + z) sinat .(13) 

xis horizontal axis,zis vertical axis where z = Oat the 
surface of still water level, /time, Gwave constant, k 
wave number, a = angular frequency, T wave period 
andh still water depth. 


The equation was formulated at flat bottom condition, 
however Hutahaean (2008) found out that the effect of 
slopping bottom on velocity potential is small, only on its 
hyperbolic term, i.e. 

Flat bottom: 

pkih+z) _|_ g-klh+z) 

coshk ih + z) =- 


Slopping bottom: 

fiiz) = ae kCh+z) + e ~ kih+z) 

Wherecris a coefficient that is a function of bottom slope 
(equation 14). It is seen thata « l.Therefore, (13) can be 
performed at sloping bottom where there will be values of 

a k ,a g 

—and—. 



1 + 


dh _dh\ 

' dx , 1 dx \ 
dh T- dh I 
1 -— 1 +—/ 
dx dx' 


(14) 


—is bottom slope. 


3.2. Wave Number Conservation Equation 
The velocity potential equation (13) is obtained from 
variable separation method, where velocity potential is 
considered as multiplication between 3 functions, i.e. 
0ix, z, t) = Xix)Ziz) Tit) , A(;c)is just an ^function 
,Z (z)is just a zfunction and T it) is just a time function. At 
(1), Ziz) = coshfcQi + z). If (13) is performed at sloping 
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bottom -=-= sinh kyh + z) -= 0 , 

3x 3x 3x 

in this equation the one with the value of zero is, 

^^=0 .(15) 

3x 

for all z value. Therefore the value of kih + z) = c, 
where cis constant, i.e. the same for the entire flowfieldof 

the wave moves. If(3) is performed on z = 0, then^- = 

3x 

Oor, 

3k k 3h , 

— =..(16) 

3x h 3x 

With (16), derivative equations higher than wave number 
can be formulated, for example for z = 0 , by ignoring 
a 2 h 
a* 2 ’ 

a 2 k 2 k /a h\ 2 

ax 2 h 2 Va xJ ' 

From this point onward, the calculation of —and—prefers 
r ax ax 2 

toz = 0. With (17) the third differential can be obtained, 
and so forth. Based on (15), the following relations apply, 

tanhkih + rf) = tanhk 0 ih 0 + 77 ) = 1 .(18a) 

coshkih + rf) = coshk 0 ih 0 + rf) .(18b) 

sinhk Qi + rf) = sinhk 0 ih 0 + rj) .(18c) 

coshk 0 ih 0 + rf) = sinhk 0 ih 0 + rf) .(18d) 

Wherer/s the water surface elevation. Therefore, based on 
(18a-d), equations containing the three elements are 
elements with values similar to the value in deep water. 

3.3. Energy Conservation Equation 

From velocity potential (1) horizontal— x velocity 

equation is obtained 

50 / 5G \ 

u = -= Gksinkx - coskx 

5x V 3 x ) 

coshk ih + z)sinot .(19) 

/ 5k 3G 

— = Gk 2 coskx + G —sinkx + 2 —ksinkx 
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5x 


5x 


5x 


a 2 G \ 

- -coskx 

5x 2 J 

coshk ih + z) sinot .( 20 ) 

and vertical—z velocity equation, 

50 

wix, z,t) = — —— = —Gkcoskxsinhk (h + z )sinot 


3z 


•( 21 ) 


— = — Gk 2 coskxcoshkih + z) sinot .(22) 

3z 

Substitute equations (20) and (22) to continuity equation 

— + — = Oand performed at the condition coskx = 

ax a z r 

sinkx = cosot = sinot = —and z =n = -, then the 

2 ' 2 

equation is divided by coshk (h + j), to obtain. 


ak a g a 2 c 

G -h 2k - - = 0.. 

3x 3x 3x 2 


..(23) 


This equation is another form of energy conservation 
equation. This equation is a relation between Gand—. The 

3x 


simplest way is by performing the assumption of a long 


wave where T^can be ignored, and in this case the 

.(24) 


•(25) 


following equation is obtained, 
ac _ g ak 

3x 2k 3x ’ 

(23) can be written as, 

a 2 c at ac 

— = G-h 2k — . 

3x 2 3x 3x 

(25) is differentiated twice against horizontal-xaxis and 

g2 q 

substituted to the term —— , and an assumption is 

3x z 

g4 q 

performed that^jis a very small number that is 

considered to be equal to zero which produce, 
ac 
ax 


\iG .(26) 

(a 3 k tak \ 2 a 2 k 2 ak) 

+2k ^c2 + 4k te) 

( a 2 k ak rr 

(4—+16k—+ 8 k 3 J 


F t a2fc ak .(^2) 

4 i^ + 16 k il + 

Therefore particle velocity equation at horizontal- 
xdirection becomes 

u = Giksinkx — \icoskx)coshk (h + z)sinot 
....(28) 

IV. DISPERSION EQUATION 

At the potential velocity equation (13), there are 2 (two) 
unknowns, i.e. energy constant Gand wave number fc; 
therefore, two equations are needed to calculate those two 
unknowns. Governing equation for analyzing the two 
unknowns are kinematic free surface boundary condition 
and momentum equation, hi its movement from the deep 
water to shallower water, evolution or transformation of 
the two unknown values will happen. The evolution is 
arranged by wave number conservation equation (15) and 
energy conservation (23) or (25). The two conservation 
equations are absorbed to the two governing equations. 

4.1 Kinematic Free Surface Boundary Condition 
Using total derivative equation(12), kinematic free 

surface boundary condition becomes w v = y~^ + u 71 ~~ ■ 
Substitute (21) , (28) and r](x,t) = Acoskxcosot and the 
equation is performed at the condition coskx = sinkx = 
cosot = sinot, 

G | k tanhk (h. + — j — yoA — (fc — p) 


coshk [h + ^ — yoA = 0 ....(29) 


4.2. Horizontal Momentum Equation 
For a function/ = fix, z, t ), where the main change is in 
the direction of horizontal-x axis, then ( 12 ) can be 
performed to obtain total acceleration equation, and 
horizontal-* and vertical-z total velocity equations are 


Du 

cLt 


3u 3u 3u 

y — +u— + w— .(30) 
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Dw 3w 3w 3 w 

— = Y — + U—+W— .(31) 

dt 'at ax az 

respectively, whereuis water particle velocity at 
horizontal-xdirection and wis particle velocity vertical-z 
direction. With (30) and (31), then Euler momentum 

equation becomes, 
au au au i a p 

Y -I - u -h w — ..(32) 

3t 3x 3z p cix 

aw aw aw i a p 

Y - 1 - u - h w — = - q .(33) 

'at ax az paz a 

At (33) the characteristics of irrational flow is performed 
at space differential, y = y, and integrated against 
vertical-z axis and dynamic free surface boundary 
condition is performed where p t] =0, pressure equation is 
obtained, i.e. 


(u 2 + w 3 ) 


p r ’ 1 3 w 1 1 

p ) 7 at 2 v 7 111 2 

+g (rj - z) 

The pressure equation is differentiated against horizontal- 
xaxis and substituted to (32) where at (32) the 
characteristics of irrational flow is performed 

3u a r v a w /i a . . a77 \ 

r m = ~ r Tx l m dz ~ li5T + < + 9 Tx) 

The completion of — fj — dzwill be done using potential 


3 r ?] 3 w 

3* z at 

velocity theory of the linear wave theory(21). 

— = -Goksinhkih + z) coskxcosot 
3 1 

77 r 

— dz = —Gcycoshkih + rj) — coshkih + z)J 
3 1 


l 


coskxcoscrt 

— [ —dz = Go(coshk(h + rf) — coshk(h + z)) 

5 x J z at 

C ksinkx — \icoskx) cosot 

From (28), 

— = Goiksinkx — | icoskx)coshk(h + z)cosot, ,hence 


a rv aw . a u v an , . 

— —dz = -, honzontal-xmomentum equation 

ax z at at at ^ 

becomes 


a u„ 


+ . (34) 


4.3. Simple Dispersion Equation 

To obtain a simple dispersion equation, convective 
acceleration at (34) is ignored, 


Y 




(35) 


Substitute (28) and water surface equation rj(x,t ) = 
Acoskxcosot, and the equation is performed at the 
condition coskx = sinkx = cosot = sinot 


YGo(k — \i)coshk (h + j) = gkA ....(36) 


Equation (29) is written as an equation for Gand 
substituted to (36), 


Y 2 o 2 (k — p) = gk yk tanhk (h. + — j — (fc — p) 

.(37) 

If the bottom slope is ignored, then (37) becomes 

Y 2 o 2 = g (k tanhk {h + - j — .(38) 

If wave amplitude is considered as a very small number, 
both to water depth and wave length, (38) becomes 

Y 2 o 2 = gktanhkh .(39) 

Then ify = lis taken, (39) becomes 
er 2 = gktanhkh .(40) 

(40) is a dispersion equation of linear wave theory (Dean, 
1991). 

Dispersion equations (37), (38), (39) and (40) have not 
met wave number conservation equation. At (37) wave 
number conservation equation (18a) is performed, hence 

Y 2 o 2 (k - p) = gk(k - (fc - p) (y)).(41) 

(41) is used to calculate wave number at the deep water. 
The dispersion equation at the shallow water is obtained 
by substituting wave number conservation (15) that can 
be stated as 

Keeping in mind that tanhk 0 ( h 0 + y) = 

lwherefc 0 ( h 0 + yj = xpn. this research used ip = 1.1, 
wheretanh (l.lir) = 0.998009, 

k(h + -) = xpTT ,ory = xpTT — kh .(42) 

Substitute (42) to (41), 

Y 2 o 2 (k — p) = gk{k — (.k — pX i/ot — kh))... (43) 

This equation is dispersion equation at the shallow water. 
However,calculation with (43) should be performed 
consecutively from deep water depth. To obtain deep 
water depth,fc 0 is calculated with (41), then deep water 

depth h 0 is the deepest between k 0 (h 0 + y) = ipn and 

— < 0.10. For water depth more than h 0 the wave 
2 ho 

number conservation equation can’t be applied. 
4.4.Complete Dispersion Equation 

hi this complete dispersion equation, the surface 
momentum equation is used completely and the wave 
number conservation equation is applied. The resulted 
equation is for calculating wave number at the shallow 
water only. Substitute (42) to (29) the first f(k,G ) = 
Oequation is obtained. 

f x (k, G) = G(k — (.k — \i)(xpn — kh)) 

coshk {h + j) — yoA = 0.(44) 

The second equation is surface momentum equation, i.e. 
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a u v /i a , an\ 

f 2 (.k,G ) = y - F (- (u 2 + w 2 ) + g —) = 

' 2 1 at \2 5x y 11 11J w a x) 


. (45) 

Where, 

Y = yGa(k — \i)coshk [h + yj ....(46a) 

Ujj- ^ = ^G 2 fc 2 (fc — | x)cosh 2 k ( h + ....(46b) 

w, ; = — ^ G 2 k 2 (.k — \i)sinh 2 k(h + -) ....(46c) 

g — = —2g(xlm — kh) ...(46d) 

At (46d) Mis substituted with wave number conservation 
equation (42). Keep in mind that based on wave number 

conservation equation, coshk (h + j)is constant number, 
i.e. coshk (h + j) = coshk 0 [h 0 + —wherefc 0 is 
calculated with (41) and deep water depth h 0 is the 
deepest between(i 0 = (ipn — ^^)^-and < 0.10. The 

values of fcandGcan be obtained by completing (35) and 
(36) with Newton-Rhaps on method, with the inputs wave 
period, wave amplitude and water depth. 


V. THE ADJUSTMENT OF y VALUE 

The value of y = 3from the previous analysis is 
theoretical value based only on Laplace equation solution. 
In this part, the adjustment of yvalue will be done using 
observation on deep water wave height. The adjustment is 
done using the relation between deep water wave height 
and wave period from Silvester (1974) and from Wiegel 
(1949 and 1964). 

By ignoring bottom slope, then (41) which is a dispersion 
equation at deep water, becomes a quadratic equation for 
k. 

y 2 a 2 = gk(\ - (“)) .(47) 

This equation has a solution if the determinant value is 
D > 0, where 

D = g 2 — 4 (y 2 <r 2 ). 


8 

2,58 

7,12 

8 

0,32 

9 

3,26 

8,02 

9 

0,32 

10 

4,03 

8,91 

10 

0,32 

11 

4,88 

9,8 

11 

0,32 

12 

5,8 

10,69 

12 

0,32 

13 

6,81 

11,58 

13 

0,32 

14 

7,9 

12,47 

14 

0,32 

15 

9,07 

13,36 

15 

0,32 


The result of the calculation on Table (3) was done using 
the value of y = 2.483 . Wavelength L on Table (3) was 

calculated using (47) obtained that Hl ” ax for all reviewed 
wave period is 0.318 or-, where it is in accordance with 

7T 

the analytical result ( Hutahaean (2019)) i.e. breaking 
occurs when^- = -, H b breaker height andL h breaker 

length. Therefore H max on column 2 is deep water wave 
height maximum for wave period on column 1, where the 
wave period is similar to 7^^ and is close enough to 
T si; that was calculated using H max on column 2. Therefore 
a conclusion can be made that the value of y = 2.483is a 
quite good value, and the maximum deep water wave 
height H 0 for wave period on column 1 is on column 2. 

VL EXAMPLE OF THE RESULT OF 
WAVELENGTH CALCULATION 

The example of the result of wavelength Lcalculation 
wave with wave period of 8 second, withv4 0 = 0.6 m 
and— = —0.01 is shown on Fig.l., Fig 2. and Fig.3. 


Fig.l shows the comparison between wavelength (40), 
(39), (38) and (43) where it is seen that (40) as dispersion 
equation of linear wave theory produces wavelength that 
is much longer than the three comparing equations. 
Wavelength (39), (38) and (43) look close, but further 
information can be seen on Fig.2. 


ForD = 0, A max = Z^With the value of H max = 
2i4 maa; , wave period is calculated from empirical 
equations of Silvester (1974), T sil = ^|l9.68//iy 3 and 

Wiegel equation (1949 and 1964), T wieg = 15.6(—j , 

H m is maximum deep water wave heighLgis the force of 
gravity. As Hi^andH m ,H max is used. 

Table.3: Wave height maximum at deep water, aty = 



oa 40 


20 


5 o 



0 2 4 6 8 10 

Water Depth h (m) 


2.483. 


T 

Hmax 

Tsu 

Twieg 

Hmax 

(sec.) 

(m) 

(sec.) 

(sec.) 

L 

6 

1,45 

5,34 

6 

0,32 

7 

1,97 

6,23 

7 

0,32 


eq. (40) eq. (39) 

eq. (38) eq. (43) 

Fig.l: Comparison between wavelength (40), (39), (38) 
and (43). 
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20 



- eq. (39) - eq. (38) 

-eq. (43) - eq. (44+45) 

Fig.2: Wavelength (39), (38), (43) and (44+45) 

Fig. 2 shows the result of the calculation using (39), (38), 
(43) and (44+45). At shallow water, wavelength from 
(39) and (38) looks much bigger than from (43) and 
(44+45). Whereas at deep water, (38) produces 
wavelength similar to that of (43) and (44+45). This 
shows that wave number conservation has a major role in 
the transformation of wavelength at shallow water, where 
at (39) and (38) wave number conservation equation is 
not performed. In addition, the changes in wavelength 
from (43) and (44+45) look linear which shows that the 
changes in wavelength as a result of water depth changes 
is dominated by wave number conservation equation (15). 
Between (39) and (38), there is a relatively big difference, 
where at (39) there is no wave amplitude as its variable as 
it is with (38). This shows that the effect of wave 
amplitude on wavelength is shortening wavelength. To 
study the effect of wave amplitude on wavelength, (38) is 
performed with different wave amplitude, i.e. 0.30 m and 
0.60 m, with the result as presented on Fig. 3, which 
shows that wavelength from a wave with wave amplitude 
0.30 m is longer that wavelength of a wave with wave 
amplitude of 0.60 m 

20 



0 

0 2 4 6 8 10 

water depth h (m) 

- wave amplitude 0.3 m - wave amplitude 0.6 m 

Fig.3: Wavelength (38) with different wave amplitude A 


VIL CONCLUSION 

This research concludes that at a space and time function, 
there is a time scale coefficient at total change or total 
acceleration. For a function f(x,z, t)with the main 
direction of change at axis — x direction, the total 
acceleration coefficient has a value of 2.483. The 
application of total acceleration equation with time scale 
coefficient at wavelength analysis produces wavelength 
that fits with the one exists in the nature. 

There are 3 factors affecting wavelength, i.e. total 
acceleration equation, wave number conservation law and 
wave amplitude. However, the main factors are the first 
and the third factors. Total acceleration plays a role in 
determining wavelength as a whole i.e at deep water and 
shallow water, wave number conservation equation plays 
a role in the transformation of wavelength at the change 
of water depth at shallow water. With the presence of 
wave amplitude effect on wavelength, the correct 
wavelength analysis is if it is performed together with 
shoaling analysis. 

Convective acceleration term at the momentum equation 
is shortening wavelength although it is relatively small. 
For practical purposes dispersion equation formulated 
without taking into account convective acceleration can 
be used. 

Wavelength research with physical model has never been 
done before. Considering that the truth of a wave theory is 
also shown by the produced wavelength, therefore the 
availability of wavelength data as the result of physical 
model is highly needed. 
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